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ABSTRACT 

An  abstract  approximation  framework  for  the  identification  of  nonlinear  distributed  parameter 
systems  is  developed.  Inverse  problems  for  nonlinear  systems  governed  by  strongly  maximal  mono¬ 
tone  operators  (satisfying  a  mild  continuous  dependence  condition  with  respect  to  the  unknown 
parameters  to  be  identified)  are  treated.  Convergence  of  Galerkin  approximations  and  the  corre¬ 
sponding  solutions  of  finite  dimensional  approximating  identification  problems  to  a  solution  of  the 
original  infinite  dimensional  identification  problem  is  demonstrated  using  the  theory  of  nonlinear 
evolution  systems  and  a  nonlinear  analog  of  the  Trotter-Kato  approximation  result  for  semigroups 
of  bounded  linear  operators.  The  nonlinear  theory  developed  here  is  shown  to  subsume  an  existing 
linear  theory  as  a  special  case.  It  is  also  shown  to  be  applicable  to  a  broad  class  of  nonlinear  elliptic 
operators  and  the  corresponding  nonlinear  parabolic  partial  differential  equations  to  which  they 
lead.  An  application  of  the  theory  to  a  quasilinear  model  for  heat  conduction  or  mass  transfer  is 
discussed. 
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1.  Introduction 


In  this  paper  we  develop  a  general  abstract  approximation  framework 
for  the  identification  of  nonlinear  distributed  parameter  evolution  systems. 
Our  intent  is  to  define  relatively  straightforward  and  easily  verified 
criteria  that  are  applicable  to  broad  classes  of  nonlinear  systems;  these 
criteria  will  guarantee  the  convergence  of  solutions  to  a  sequence  of  finite 
dimensional  Galerkin  approximation  based  parameter  estimation  problems  to 
a  solution  of  the  original,  underlying,  infinite  dimensional  identification 
problem.  The  results  that  we  present  below  generalize  and  extend  the 
theory  recently  developed  by  Banks  and  Ito  in  [2]  and  [3]  for  regularly 
dissipative  or  abstract  parabolic,  linear  systems.  It  is,  to  the  best  of  our 
knowledge,  the  first  such  general  approximation  theory  for  inverse 
problems  involving  nonlinear  distributed  systems. 

The  sufficient  conditions  set  down  in  our  framework  include  a 
relatively  mild  continuity  assumption  with  respect  to  the  unknown 
parameters  to  be  identified,  an  equi-boundedness  and  an  equi-strong 
monotonicity  assumption  on  the  nonlinear  operator  describing  the  system 
dynamics.  In  addition  our  theory  requires  a  standard  approximation 
assumption  on  the  Galerkin  subspaces  used  to  effect  the  finite  dimensional, 
or  finite  clement,  approximations.  We  demonstrate  that  solutions  to  the 
finite  dimensional  identification  problems  approximate  a  solution  to  the 


infinite  dimensional  identification  problem  via  a  convergence  result  for 


solutions  to  the  forward  problems, 
of  nonlinear  evolution  systems  and 


This  result  is  obtained  using  the  theory  i  For 


a  nonlinear  analog  of  the  well-known 


Trotter-  Kato  approximation  result  for  linear  semigroups. 
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In  the  present  paper,  we  are  concerned  only  with  theory; 
implementation  questions  and  conclusions  drawn  from  our  numerical  or 
computational  studies  will  be  reported  on  elsewhere.  Also,  while  we  have 
tried  to  make  our  framework  as  versitile  as  possible,  the  treatment  below 
does  have  limitations.  For  example,  our  theory  can  handle  quasi-autonomous 
systems  but  it  is  not  applicable  in  the  fully  nonautonomous  case.  The 
development  of  a  general  theory  which  can  handle  nonlinear  systems 

involving  time  dependent  operators  requires  additional  effort  and  is 

currently  the  focus  of  our  ongoing  investigations.  The  particular 
difficulties  inherent  in  the  time  dependent  case  will  be  described  in 
greater  detail  in  our  discussions  below. 

We  provide  a  brief  outline  of  the  remainder  of  the  paper.  In  Section 
2  we  state  a  fundamental  existence  and  uniqueness  result  for  infinite 

dimensional  nonlinear  systems  and  prove  a  general  approximation  result 
which  is  especially  well  suited  for  application  in  the  context  of  the  inverse 
problems  which  are  the  central  focus  of  our  study.  In  Section  3  we 
define  a  class  of  nonlinear  distributed  systems  and  the  associated 
parameter  identification  problems.  We  define  the  Galerkin  approximations 
and  prove  the  general  convergence  result.  Section  4  contains  some  examples. 
Wc  show  that  our  nonlinear  theory  subsumes  the  linear  theory  presented  in 
[2]  and  [3]  as  a  special  case;  we  also  consider  the  application  of  our 

framework  to  a  class  of  nonlinear  elliptic  operators  and  the  corresponding 
nonlinear  parabolic  partial  differential  equations  to  which  they  lead.  In 
particular,  we  look  at  the  application  of  our  results  to  a  well  known 
quasilincar  model  for  heat  conduction  or  mass  transfer.  In  Section  5  we 
summarize  our  findings  and  provide  some  concluding  remarks. 


2.  An  Approximation  Result  for  Nonlinear  Evolution  Systems 


Let  Xc  be  a  Banach  space  with  norm  |  |0.  We  consider  the  nonlinear, 
quasiautonomous  initial  value  problem  in  XQ  given  by 

(2.1)  x0(t)  +  A0x0(t)  B  f0(t),  0  <  t  S  T, 

(2.2)  x  Q(0)  =  x° 

where  x°  €  XQ,  fQ:  [0,T]  -  XQ  and  the  nonlinear  operator  AQ:  XQ  -  2  0  is  in 
general  multivalued,  not  everywhere  defined,  and  not  continuous.  The 
existence  of  solutions  to  the  initial  value  problem  (2.1),  (2.2)  and  the 
subsequent  approximation  result  to  follow,  are  both  consequences  of  Theorem 
2.1  to  be  given  below. 

We  shall  require  the  following  definitions.  Let  X  be  a  Banach  space  with 
norm  |-|x.  For  A:  X  -  2X,  a  nonlinear,  multivalued  operator,  the  domain  and 

range  of  A  are  defined  by  Dom(A)  =  {x  €  X:  Ax  *  0}  and  R(A)  =  u  Ax 

XCdom(A) 

respectively.  We  say  that  the  operator  A  is  accretive  if  for  every  X  >  0, 
xrx2  €  Dom(A)  and  yt  €  Axp  y2  €  Ax2  we  have 

lxi  -  x2jx  «  |xt  -  x2  +  X(yt  -  y2)|x. 

We  say  that  A  is  m-accretive  if  A  is  accretive  and  R(I  +  XA)  =  X  for  some 

X  >  0.  We  note  that  if  A  is  m-accretive  then  R(I  +  XA)  =  X  for  every  X  >  0 

and  for  each  X  >  0  the  resolvent  of  A  at  X,  J(X;A):  X  -*  X,  a  single  valued, 
everywhere  defined,  nonlinear  operator  on  X  can  be  defined  as  J(X;A)  = 
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A  two  parameter  family  of  nonlinear  operators  (U(t,s):  0  i  s  i  t  <  T) 
defined  on  a  subset  0  c  X  is  called  a  nonlinear  evolution  system  on  fi  if  for 
each  x  G  0  we  have  U(t,s)x  G  A,  U(s,s)x  =  x  and  U(t,r)U(r,s)x  =  U(t,s)x  for 
0£sSr«t«T  and  U(t,s)x  is  continuous  from  the  triangle  A  =  {[s,t]: 
0  ^  s  S  t  S  T)  into  X. 

A  strongly  continuous  function  x:  [0,t]  -  X  is  called  a  strong  solution  to 
the  quasiautonon.ous  initial  value  problem 


(2.3)  x(t)  +  Ax(t)  3  f(t),  0  <  t  S  T 


(2.4)  x(0)  =  x° 


where  f:  [0,T]  -  X  and  x°  €  X  if  x  is  absolutely  continuous  on  compact 
subintervals  of  (0,T),  differentiable  almost  everywhere  and  satisfies  f(t)  - 
x(t)  G  Ax(t)  for  almost  every  t  G  [0,T]  and  x(0)  =  x°. 


Theorem  2.1.  Let  X  be  a  Banach  space  with  norm  ||x  and  suppose  that  A:  X  - 
2X  and  f:  [0,T]  -*  X  appearing  in  (2.3)  satisfy 

( 1 )  there  exists  an  w  G  R  for  which  the  operator  A  +  U)I  is  m -accretive, 

(2)  f  G  L1(0,T;X). 

Then  a  unique ,  nonlinear  evolution  system  (U(t,s):  0  S  s  5  t  S  T)  on  Dom(A)  can 
be  constructed  which  satisfies 


(i)  |U(t,s)4>  -  U(t,s)0|x  «  eu(t")|(J>  -  <HX,  for  p,p  €  Dom(A)  and 
0  <  s  <  t  (  T, 

(ii)  |U(s+t,s)0  -  U(r+t,r)0|x  S  2f  eu(t-T)|f(T+s)  -  f(T+r)|xdr, 

_ 

for  all  ip  G  Dom(A)  and  all  t  >  0  such  that  s+t,  r+t«T. 
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(iii)  if  x°  6  Dom(  A)  and  the  initial  value  problem  (2.3),  (2.4)  has  a  stiong 
solution  x,  then 

x(t)  =  U(t,s)x(s),  for  0  (  s  U  5  T. 

When  x°  €  Dom(  A),  the  strongly  continuous  function  x:  [0,T]  -*  X  given  by  x(t) 
U(t,0)x°  is  referred  to  as  a  mild  or  generalized  solution  to  (2.3),  (2.4). 

Theorem  2.1  is  a  direct  consequence  of  results  given  by  Crandall  and 
Evans  and  Evans  in  [7]  and  [9].  Henceforth,  we  shall  assume  that  AQ:  XQ  -  2xo 
and  fQ:  [0,T]  -  X,  satisfy  (1)  and  (2)  in  the  statement  of  Theorem  2.1  and 

that  xQ  e  Dom(A0).  We  then  let  (U0(t,s):  0  S  s  ^  t  «  T)  denote  the 
corresponding  nonlinear  evolution  system  on  Dom(A0)  and  consider  the 
approximation  of  mild  solutions  to  the  initial  value  problem  (2.1),  (2.2). 

Our  approximation  result  is  in  the  spirit  of  those  given  for  nonlinear 
semigroups  and  evolution  systems  by  Crandall  and  Pazy  in  [8]  and  Goldstein 
in  [10],  However,  our  theorem  differs  from  these  earlier  treatments  in  two 
ways.  First,  we  require  that  the  time  dependent  perturbation  fQ  be  only  as 
opposed  to  it  being  continuous  as  in  [8]  and  it  satisfying  a  Lipschitz-like 
condition  in  [10].  This  distinction  is  especially  relevant  in  the  case  of  control 
systems  where  discontinuous  input  is  common.  The  second  difference  is  that 
we  give  our  result  in  a  form  that  is  most  appropriate  for  application  to  the 
development  of  a  general  approximation  theory  or  framework  and 
computational  schemes  for  the  parameter  identification  problems  to  be 
discussed  in  the  next  section. 

We  shall  require  some  set  theoretic  notation.  For  sets  Hn,  n=0,l,2,...,  by 
lim  HnDH0  we  shall  mean:  Given  xQ  €  HQ,  there  exist  xn€ Hn  such  that  xn-xQas 


Theorem  2.2.  For  each  n  e  Z+  =  {1,2,3,  ...}  let  Xn  be  a  closed  linear  subspace  of  XQ. 

x 

For  n  =  0,1 . let  An:  Xn  -  2  n  be  a  possibly  multivalued  nonlinear  operator  on  Xn,and 

let  f  :  [0,T  ]-  Xn  bean  Xn-valued  measurable  function  defined  on  [0,T],  Suppose 
that  there  exists  an  uQ  e  R,  independent  of  n,  for  which  the  operators  An  +  u0I 
are  m -accretive,  that  there  existsa  function  g€L1(0,T;Xo)  for  which  |f n ( t )|  ?  g(t), 
a.e.te  [0  ,T],  and  that  lim  Dn  3  DQ  where  Dn  =  Dom(An)  and  DQ  =  Dom(A0). 
Suppose  further  that  for  some  XQ>0we  have 

(2.5)  lim  J(X0;  An  +  uQI)0n  =  J(X0;  AQ  +  uQI )4>0 
whenever  <t>  €  X  with  lim  <$>  =  $n€  X., 

n  n  *oo  n  u  u 

and  that 

lim  f  (t)  =  f_  (t)  for  a.e.te  [0  ,T]. 

n-°°  n  u 

Then  for  each  n  €  Z+  there  exists  a  unique  nonlinear  evolution  system  (Un(t,s): 
0  S  s  U  i  T)  on  corresponding  (in  the  sense  of  Theorem  2.1)  to  An  and  fn 
and  for  <t>  €  D  with  lim  <t>  =  <t>n  €  Dn  we  have 

n  n  n-*00  n  0  u 

(2.6)  Hm  Un(t,s)0n  =  U0(t,s)4>0,  0  S  s  «  t  S  T, 

with  the  limit  being  uniform  in  t  for  t  e  [s,T], 

Proof.  We  follow  Goldstein  (see  [10],  [11])  and  use  an  approach  first  suggested 
bv  Kisynski  [13]  for  demonstrating  the  convergence  of  approximations  to 
linear  semigroups,  to  prove  the  theorem  via  an  application  of  our  existence 


result.  Theorem  2.1. 


Let  X  =  {x  =  {x  }„_n:  x  G  X  ,  n  =  0,1,2 .  and  lim  x  =  xn}  and  for 

n  ’  fi  — u  n  n  *oo  n  u 

A  A 

x  G  x  set  llxll  =  sup|xn|Q.  Then  II  •  II  defines  a  norm  on  the  linear  vector  space 

n 

X,  and  the  space  X  together  with  the  norm  II  S  is  a  Banach  space.  Define  the 
operator  A:  x  -  2X  by 

dom(^)  =  {x  =  {xn}“=0  ex:  xn  e  Dom(An),  and  for  each  n  =  1,2,  ... 
there  exists  ay  e  A  x  such  that  lim  y  =  yA  e  A^x.}, 

J  n  n  n  n-*00  n  ^  0  0 

for  x  e  Dom(/l),  y  =  {yn}~=0  e  /lx  if  and  only  if  yn  €  Anxn  , 
n  =  0,1,2,  ...  and  lim  y  =  yn. 

n-“  '  n  U 


Define  an  essentially  x-valued  function  f  on  the  interval  [0,T]  by  f(t)  = 
{fn(t)}”_0.  The  assumptions  on  the  f  are  such  that  fn(t)  -  fQ(t)  for  almost 
every  t  €  [0,T].  However,  by  appropriately  redefining  on  a  set  of  measure 
zero,  we  may  infer  from  the  assumptions  on  the  functions  f  that  f:[0,T]  -  x 
with  f  G  LjfO.TiX). 

It  is  readily  seen  that  the  operator  A  +  uQI  is  m-accretive.  Let  x1  = 

<xX=0'  *2  =  {x^>r=o  e  Dom(zl)  and  ict  y1  =  {y*}”=o  e  A and  V2  =  <y’C=o  e 
Ax  .  Since  for  each  n  =  0,1,2,  ...,  An  +  u)QI  is  assumed  to  be  m-accretive,  for 
X  >  0  we  have 


lx1  -  x2 H  =  sup  |xx  -  x2L  «  suplx1  -  x2  +  X(y*  +  ui-x1  -  (y2  +  wr,x2))|„ 

^  '  n  n'O  n  n  w  n  On  wn  0  n7/,0 

n  n 

=  Ik1  -£2  +  X(y 1  +  u)^1  -(y2  +  u0x2))ll. 


•  •  ^  CD  A 

and  therefore  that  A  +  wQI  is  accretive.  Now  let  y  =  {yn}n=0  e  x  and  set  x  = 
{xn}”_Qwith  xn  =  J(XQ;  An  +  uQI)yn,  n  =  0,1,2,  ...  where  XQ  is  chosen  as  in  (2.5).  It 
is  immediately  clear  that  for  each  n  =  0,1,2,  ...,  x  e  Dom(A  )  C  X  .  Since 


••M 


'<  > 


€  x  we  have  lim  yn  -=  y0  and  therefore,  by  assumption  (2.5),  that  lim  xn  =  x0 
or  x  G  x.  Setting  zn  =  (yn  -  (1  +  X0ta0)xn)/X0,  n  =  0,1,2,  ...,  it  follows  that 
zn  €  Anxn  and  lim  zn  =  zQ  G  A0xQ.  We  conclude  that  x  G  Dom(/l),  (I  + 
X0(^+w0I))x  3  y,  and  that  R(l  +  XQ(/1  +  uQl))  =  x. 

We  have  shown  that  the  operator  A  and  the  function  f  satisfy  conditions 
(1)  and  (2)  given  in  the  statement  of  Theorem  2.1.  Therefore,  a  unique 
nonlinear  evolution  system  {{/(t,s):  0  S  s  t  t  t  T)  on  Dom(/l)  corresponding  to 
A  and  f  can  be  constructed  with  (7(t,s)  =  {Un( t,s))”_0.  Using  assumption  (2.5) 
it  can  be  shown  that  Dom(/l)  =  (x  =  {x  .  G  x:  x  6  D  .  n  =  0,1,2,  ...  and 
lim  xn  =  xQ).  Since  £(t/(t,s))  C  x,  it  follows  that 

(2.7)  lim  C/n(t,s)0  =  f/Q(t,s)0  0  S  s  S  t  «  T 

n-*00 

w'hcncver  0  e  D  and  lim  <5  =  <pn  €  D..  Since  each  of  the  operators  A  and 

n  n  ^  — *co  n  u  u  n 

the  functions  fn  satisfy  conditions  (1)  and  (2)  of  Theorem  2.1,  unique 
nonlinear  evolution  systems  (Un(t,s):  0  5s  U  t  T)  on  Dn  corresponding  to  An 

_  CD  _ 

and  f  can  be  constructed.  Recalling  that  Dom(^)  C  x  D  ,  we  may 

"  n=0  n 

define  the  family  of  operators  (V(t,s):  0  5  s  5  t  S  T)  on  Dom(  A)  by 

(2.8)  V(t,s)x  =  {Vn(t,s)xn}"=0  =  (Un(t,s)xn}”=0 

for  x  =  {xn}"=0  G  Dom(^).  Uniqueness  (see  [9])  dictates  that  for  each  n  = 
0.1.2,  ...,  f/n(t,s)xn  =  V Jt  s)xn  whenever  (xn}“=0  €  Dom(/t).  This  together 
with  (2.7)  and  (2.8)  establish  (2.6).  The  fact  that  the  convergence  in  (2.6) 
is  uniform  in  t  for  t  €  [s,TJ  is  argued  exactly  as  it  was  for  the 
convergence  of  approximations  to  nonlinear  semigroups  in  the  proof  of 
Theorem  3.2  in  [10]. 


3.  An  Approximation  Theory  for  Identification  Problems 


Let  H  be  a  real  Hilbert  space  with  inner  product  <-,•>  and 

corresponding  norm  |  |.  Let  V  be  a  reflexive  real  Banach  space  with 
norm  I!  •  II  and  let  V*  be  its  dual.  (All  of  our  theory  can  be  developed  in 
complex  spaces  if  necessary;  see  (6).)  We  denote  the  usual  dual  norm  on  V* 
by  II  II.  and  assume  that  V  is  densely  and  continuously  embedded  in  H 
with  |v|  S  trllvll,  v  e  V,  for  some  positive  constant  jc  Identifying  H  with 

its  dual,  we  obtain  V  C  H  =  H*  C  V*.  For  0  6  V*  and  v  e  V  the  duality 

pairing  between  0  and  v  is  denoted  by  <0,v  >.  When  0  €  H,  its  pairing 

with  v  €  V  agrees  with  the  inner  product  of  0  with  v.  It  follows  for  u  e 
H  and  v  e  V  that  Hull.  S  m|u|  and  Bvll.  $  <z2llvll.  Let  Q  and  Z  be  metric 
spaces  and  let  Q  be  a  nonempty,  sequentially  compact,  subset  of  Q.  The 
spaces  Q  and  Z,  and  the  set  Q  are  referred  to  as  the  parameter  space,  the 
observation  space,  and  the  admissible  parameter  set  respectively. 

We  recall  that  a  single  valued  operator  A:V-.V*  is  hemicontinuous  if 

lim  A(u  +  tv)  =  Au  for  all  u,veV  where  the  limit  is  taken  in  the  weak  sense. 

t-o 

For  each  q  €  Q  let  A(q):  V  -  V*  be  a  single  valued,  hemicontinuous,  (in 
general,  nonlinear)  operator  satisfying: 

(A)  (Continuity):  For  each  v  e  V,  the  map  q  -•  A(q)v  is  continuous  from 
Q  C  Q  into  V*. 

(B)  (Equi  V-  monotonicitv):  There  exist  an  u  G  R  and  an  a  >  0,  both 
independent  of  q  e  Q,  such  that 

<A(q)u  -  A(q)v,  u~v>  +-  u^u—vl2  £  allu-vll2. 


for  every  u,v  €  V. 
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Thcorcm  3.1.  For  each  q  €  Q  the  operator  AQ(q):  Dom(A0(q))  C  H  -  H  is 
densely  defined  and  the  operator  AQ(q)  +  ul  is  m -accretive. 

Proof.  We  first  show  that  for  each  q  €  Q  the  operator  A(q)  +  ul:  V  -  V*  is 

coercive.  If  {vn}  c  V  with  lint  llvnS  =  <*>,  then  from  assumptions  (B)  and  (C)  we 
obtain 

lim  <(A(q)  +  ud)v  ,v  >/ 1! v  II 

n  —*co  n  n  n 

=  lim  {( <A(q)v  -  A(q)0,v  >  +  to|v  |2)/ 1 v  II  +  <A(q)0,v  >/lv  »} 

>  H^rn  { cx II  v n  II 2/  II  v „  II  -  !<A(Q)e,vn  >|/ *  vn  I } 

?  lim  {a II v  II  -  II A(q) 0 H .}  J  lim  aJv  I  -  0  =  <*> 

n-*00  n  n~“°  n 

where  0  denotes  the  zero  vector  in  V.  It  follows  that  for  each  V  >  0,  the 
operator  I  +  X(A(q)  +  ul):  V  •*  V*  is  monotone,  everywhere  defined  on  V, 
hemicontinuous,  and  coercive.  Consequently  R(I  +  X(A(q)  +  ul))  =  V*  (see 
Barbu  (6],  Theorem  II.  1.3)  and  therefore  X(I  +  x(A0(q)  +  ul))  =  H.  Also,  for 
u,v  €  Dom(A0(q)),  we  may  use  assumption  (B)  to  conclude 


[l  +  ]tu-v|2  S  |u~v|2  +  Xoell 


<  | u  v | 2  +  X<(A(q)  +  ul)u  -  (A(q)  +  ul)v,  u-v> 

=  <(I  +  X(A(q)  +  ul))u  -  (I  +  X(A(q)  +  wl))v,u-v> 
S  |(I  +  X(AQ(q)  +  wl)u  -(1  +  X(AQ(q)  +  ul)v|  |u-v| 


|u-v|  «  |u  -  v  +  X((A0(q)  +  ul)u  -  (A0(q)  +  wl)v)| 

which  proves  that  AQ(q)  +  ul  is  m-accretive  on  Dom(A0(q))  C  H. 

To  show  Dom(  AQ(  q) )  =  H,  we  let  u  €  H  and  for  each  n  =  1,2, ...  we  set 


un  =  J(  1 /n;AQ(q)+wl)u  €  Dom(A0(q)).  Then,  arguing  as  we  have  above,  we 
find 

|un|2  +  ( 1  / n)otP un II 2  S  <u  -  (l/n)A(q)9,un> 

«  |u|  |un|  +  (l/n)lA(q)0IJunl 

where  9  is  again  the  zero  vector  in  V.  But  then 

(3.3)  (l/2)|un|2  +  (l/n)(a/2)llunll2  «  (l/2)|u|2  +  (l/n)(l/2«)l A(q)0B2 

S  (l/2)|u|2  +  (l/n)(B2/2a), 

from  which  it  immediately  follows  that  the  un  are  uniformly  bounded  in  H. 
Indeed,  from  (3.3)  we  see  that  (l/n)Iunl2  and,  hence  lu^/vTn,  is  bounded  so 
that  11  u  II / n  —  0  as  n-** 

Also,  assumption  (C)  yields 

II u n  -  ufl,  =  (1/n) ll( A0(q)  +  wl)unl.  <  (l/n){(B  +  UM2)lunll  +  B). 

Since  the  last  term  in  the  estimate  above  tends  to  zero  as  n  -*  °°,  we  find 
un  -  u  in  V*  as  n  -*  ®.  This,  together  with  the  fact  that  V  is  dense  in  H 
imply  that  un  -*  u  weakly  in  H  as  n  -  “  from  which  Dom(AQ(q))=H 
immediately  follows. 

In  light  of  Theorem  3.1,  we  may  apply  Theorem  2.1  with  X  =  H,  A  = 
AQ(q)  and  f  =  f(  ;q).  We  conclude  that  there  exist  a  unique  nonlinear 
evolution  system  (U0(t,s;q):  OSsStSTJonH  satisfying  (i),  (ii)  and  (iii). 


The  mild  solution  uQ(  ;q):  [0,T]  -•H  to  the  initial  value  problem  (3.1),  (3.2)  is  given 
by  uQ(t;q)  =  Uo(t,0;q)u°(q)  for  t  €[0,T], 


Remark.  Under  additional  hypotheses  on  f(  - ;q)  and  u°(q)  other  existence 
results  can  be  applied  to  obtain  somewhat  different  notions  of  a  solution  to 
the  initial  value  problem  (3.1),  (3.2).  For  example  (see  [6,  p.  1 40- 1 44])  if 
f(-;q)  €  W1,1(0,T;H)  and  u°(q)  €  Dom(AQ(q)),  then  there  exists  a  unique  u(-,q): 
[0,T]  -  V  satisfying  u(  ;q)  €  W^O.TjH),  A(q)u(;q)  e  LJ0,T;H)  and  u(t;q)  + 
A(q)u(t;q)  =  f(t;q)  a.e.  t  G  [0,T],  Or,  if  u°(q)  €  H  and  f(;q)  €  L2(0,T;V*) 
then  there  exists  a  unique  u(  ;q)  which  is  V*-valued  absolutely  continuous 
almost  everywhere  on  [0,T],  u(  ;q)  €  C(0,T;H)  n  L2(0,T;V),  u(-;q)  e  L2(0,T;V*) 
and  u(t;q)  +  A(q)u(t;q)  =  f(t;q),  a.e.  t  €  [0,T].  If,  in  addition,  the  mapping  t  -♦ 
t7f'(t;q)  is  an  element  in  L2(0,T;V*)  for  some  y  >  1,  then  the  mapping  t  -* 
t^u(t;q)  is  in  L2(0,T;V)  n  Lm(0,T;H).  In  particular,  when  f(  ;q)  =  0,  the 
nonlinear  semigroup  (SQ(t;q):  0  U  t  T)  on  H  defined  by  SQ(t;q)  =  UQ(t;0;q), 
t  €  [0,T],  with  generator  ~A0(q)  behaves  like  a  holomorphic  linear  semigroup 
in  that  it  smooths.  That  is,  SQ(t;q)u0(q)  e  Dom(A0(q)),  t  €  (0,T],  and  the 
mapping  t  -*  t  ^  S(t;q)u°(q)  is  an  element  in  L,„(0,T;H)  for  every  u°(q)  €  H. 
Also,  some  generalizations  are  possible.  For  example,  in  assumption  (B).  the 
term  allu  -  v  II 2  can  be  replaced  by  a  term  of  the  form  oc( P u  — v  II ) II u  -  vll  where 
<*(•)  is  a  continuous,  strictly  increasing  function  on 
[0,")  satisfying  a(0)  =  0  and  lim  a(x)  =  ®.  Or,  the  terms  Uu  —  vil2  in  (B) 

x-a> 

and  II v II  in  (C)  can  be  replaced  by  II u  -  vllp  and  B v P p_1,  respectively,  for 
any  p  £  2. 


The  development  of  computational  methods  for  the  solution  of  the 
infinite  dimensional  optimization  problem  (ID)  requires  the  finite  dimensional 
approximation  of  the  abstract  initial  value  problem  (3.1),  (3.2).  The  general 
framework  that  we  are  proposing  is  based  upon  a  classical  Galerkin  approach. 
For  each  n  =  1,2,  ...  let  Hn  denote  a  finite  dimensional  subspace  of  H  which  is 
a  subset  of  V.  Let  Pn:  H  -*  Hn  denote  the  orthogonal  projection  of  H  onto  H 
with  respect  to  the  <-,•>  inner  product.  We  assume  that  the  approximating 
subspaces  Hn,  and  the  projections  Pn  satisfy 

(D)  For  each  v  €  V,  limSP  v  —  vlt  =  0. 

n-“  n 

Note  that  assumption  (D)  and  V  densely  and  continuously  embedded  in  H 
imply  that  lim  |P  u  -  u|  =  0  for  each  u  6  H. 

n-®  n 

For  each  q  €  Q  and  n  =  1,2,  ...  we  define  the  single-valued  operator  An(q): 
H  -  H  by  A  (q)u  =  v  for  u„  €  H  where  v  satisfies 

u  ii  n  n  n  n  n 

<A(q)un,Wn>  =  <vn,wn>,  wn  6H, 

That  An(q)  is  a  well  defined  operator  from  into  Hn  follows  from  the  Riesz 
Representation  Theorem  applied  to  the  Hilbert  space  H_  and  the  bounded 
linear  functional  <A(q)un, ->  on  Hn.  Also,  define  f n(  - ;q):  [0,T]  -  Hn  and 
u°(q)  €  Hn  by  fn(t;q)  =  Pnf(t;q),  0  5  M  T,  and  u°(q)  =  Pnu°(q),  respectively. 
Note  that  fn(  -;q)  e  L,(0,T;Hn)  c  L^O.TjH)  and  that  |fn(t;q)|  «|f(t;q)|  for  q  eQ 
and  almost  every  t€  [0,T], 

We  consider  the  sequence  of  approximating  identification  problems  given 
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(lDn)  Given  observations  z  €  Z,  determine  parameters  qn  €  Q  which 
minimize 
4>n(q)  =  <Kun(q);z) 

where  un(q)  =  un(  ;q)  is  a  mild  solution  to  the  initial  value  problem 
in  H 

n 

(3.4)  un(t)  +  An(q)un(t)  =  fn(t;q),  0  <  t  S  T 

(3.5)  un(0)  =  u°(q) 
corresponding  to  q  e  Q. 


From  the  definition  of  the  An(q)  and  the  assumptions  (B)  and  (C)  on  A(q), 
using  arguments  analogous  to  those  used  to  prove  Theorem  3.1,  it  can  be 
shown  that  the  operators  An(q)  +  ul  are  m-accretive  on  H  .  It  then  follows 
from  Theorem  2.1  that  for  each  n  =  1,2,  ...  there  exists  a  unique  nonlinear 
evolution  system  (Un(t,s;q):  0  $  s  U  t  T)  on  Hn  satisfying  (i)  -  (iii)  in  the 
statement  of  that  theorem  with  X  =  Hn,  f(t)  =  fn(t;q),  and  x°  =  u°(q).  The 
mild  solution  to  the  initial  value  problem  (3.4),  (3.5)  is  given  by  un(t;q)  = 
Un(t,0;q)u°(q),  t  €  [0,T]. 

If  we  assume  for  the  moment  that  the  approximating  identification 


exist  a  subsequence  {qn  }  C  {qn}  and  a  q  €  Q  such  that  lim  qn  =  q.  For  any 
q  €  Q  the  continuity  of  4>  implies 

0{q)  =  <t>(u0(q);z)  =  <J>  ( lim  u  (q  );z) 

j-“  nj  nJ 

=  lim  *(un  (q  );z)  =  lim  <p  (q  ) 
j-*“  nj  nj  j-40  nj  nj 

«  lim  <t>  (q)  =  lim  4>(u  (q);z) 

j-»CO  “j  j  -*CO  nj 

=  *  (lim  u  (q);z)  =  4>(u0(q);z) 

j  — *oo  "j  W 

=  <Wq). 

Note  that  in  the  discussion  above  we  did  not  assume  that  a  solution  to 
problem  (ID)  exists.  But  rather  we  have  shown  that  the  existence  of  solutions 
qn  to  the  approximating  problems  (IDn)  and  (3.6)  imply  the  existence  of  a 
solution  q  to  problem  (ID).  When  the  solution  to  problem  (ID)  is  unique,  the 
sequence  (qn)  itself  converges  to  q. 

The  existence  of  a  solution  qn  to  problem  (IDn)  for  each  n  =  1,2,  ...  will 
follow  from  the  compactness  of  Q  and  the  continuity  of  $  once  the  continuous 
dependence  result:  lim  u  fq  )  =  u  (qn)  in  C(0,T;H  )  whenever  {q  I  C  Q  with 

ni  trK>  •*  ni  n  u  n  m 

lim  q  =  q_,  has  been  established.  Although  continuous  dependence  for  the 

m-°°  m  0 

finite  dimensional  systems  (3.4),  (3.5)  could  be  demonstrated  via  a 
modification  to  any  one  of  a  number  of  familiar  continuous  dependence 
results  for  ordinary  differential  equations  (see,  for  example.  Hale  [12], 
Theorem  1.3.4),  it  is  also  easily  handled  with  the  approximation  theory 
developed  in  the  previous  section.  This  and  the  convergence  in  (3.6)  are 
addressed  in  the  following  theorem. 


Theorem  3.2.  If  assumptions  (A)  -  (D)  hold,  then 


(a)  If  {qn}  C  Q  with  lint  qn  =  qQ  then  lim  un(qn)  =  uQ(q0)  in  C(0,T;H),  and 

(b)  If  (qm)  C  Q  with  lim  qm  =  qQ  then  lim  un(qm)  =  un(q0)  for  each  n  €  Z+. 

m-*®  m-400 

Proof.  Assumption  (D)  and  the  continuity  of  the  map  q  -*  u°(q)  from  Q  into  H 
imply  lim  u°(q  )  =  u°(qn)  in  H.  Hence,  we  will  have  verified  (a)  if  we  can 

n— *oo  n  n  u 

show  that  lim  Un(t,s;qn)wn  =  UQ(t,s;q0)w0,  0  i  s  it  5  T,  uniformly  in  t  for 
t  €  [s,T]  whenever  w  e  H  with  lim  w  =  wQ  €  H.  We  argue  this  using 

n~*°° 

Theorem  2.2.  Note  that  assumption  (D)  implies  lim  H  D  H  and  assumption 

n-°°  n 

(D)  together  with  the  assumed  continuity  of  the  map  q  -*  f(t;q)  from  Qc Q 
into  H  for  almost  every  t  €[0  ,T]  imply  lim  f  (t;qj  =  f(t  ;q„)  in  H  for  almost 

n-*oo  n  n  u 

every  t€  [0,T]  with  the  fn(  ;dn)  dominated  by  a  function  g€Lx(0,T;H) 
which  is  independent  of  n.  Thus,  we  need  only  to  demonstrate 
that  for  some  XQ  >  0  we  have 

(3.7)  limJ(X0;An(qn)  +  wl)wn  =  J(XQ;A0(q0)  +  ul)w0 

n*"400 

in  H  whenever  w  €  H  ,  n  €  Z+  with  lim  w  =  w_. 

n  n  n~4<x>  n  ^ 

Let  XQ  >  0  and  set  vn  =  J(X0;An(qn)  +  ul)wn  and  vQ  =  J(XQ;A0(q)  +  wl)wQ. 
We  first  show  that  fl v^ H  is  uniformly  bounded  in  n.  From  assumption  (B)  we 
obtain 

X0allVn“2  S  +  X0<A^n)Vn  ”  A^n)0-Vn> 

=  <0  +  xo(A„(qn)^i))v„,v„>-|v/ 

+  ^o<A(qO)0  “A(Cl„)0>Vn>"  XO<A(qO)0’Vn> 

=  <Wn’Vn>  "  lVJ2  +  X0<A(qo)6  “  “  XO<A(qO)0’Vn  > 

s  lwnll.lvnl  +  XolA(qo)0  -  A(qn)0«Jvnl  +  XQ» A(qo)0 l.lvj 


9- 


where  0  denotes  the  zero  vector  in  V.  This  estimate  together  with  assumption 


(C)  yields 


l|vni  «  (xoa)"1^wnl  +  «  1|lA(qn)6  “  A(q0)9S.  +  a‘1S. 


Recalling  assumption  (A)  and  that  lim  w  =  wn  in  H,  we  find  that  the  desired 


uniform  bound  on  II v  II  has  been  established. 


Once  again,  from  assumption  (B),  we  find 


X0«||v„-v0112  «  +  X0<A(qn>Vn  “  A(<OV0’V„  “  V 

=  X0WiVn~V0l2  +  X0<A(qn)Vn  “  A(q0)VV„  “  Pn 

+  xo<A(qn)vn  -  A(q0)v0,Pnv0  -  v0> 

+  xo<A(q0)vo  ~  A(qn>vo’vn _  V0> 


=  X  u<P  v-  -  v.,v  ~  v_>  +  <(I  +  X.(A  (q  )  +  ul))v 

0  n  0  0’  n  0  v  0V  nVMn7  //  n 


(1  +  X0(A0(q0)  +  u>I))v0,vn  -  PnvQ>  +  <vQ -vn,vn-Pnv0> 


+  X0<A(qn)Vn  “  A(q0)V0’PnV0  -  V 

+  X0<A(q0)v0  -  A(qn)v0,vn  -  vQ> 

=  X0U<PnV0  "  V0,Vn  ~  V0>  +  <Wn  "  W0’V„  "  P„V  "I  Vn  ~PnV0  I* 

+  Xo<A(q„)Vn  "  A(q0)V0’PnV0  "  V0> 

+  xo<A(q0)vo _  A(q„)vvn "  V 


5  X0UUPnV0“V0ll*l,Vn“V011  +  1 1  *  *' VIT'V 1 


+  »wn-w0BJIPnv0-V()l  +  X0HA(qn)vn  -  A(q0)v0UPnv0  -  v0# 


+  x0”A(q0)v0  -  A(qn)v0ll.llvn  -  vQll. 


1  ,  n  , 

The  estimate  ab  S  —  a  +  -  b  for  any  n  >  0  and  assumption  (C)  allow  us  to 
2n  2 


argue 


^llv  -vjl2  «  —  Xq  IIP  v  -vjl2  +  II w  — w„8?  +  1  w  -wJ.lP  v„-v„ 
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+  V0UA(qn)vn  -  A(q0)v0«.IPnv0  -  vQ»  +  ^llA(q0)v0  -  A(qn)v0I? 


<  r  3u^  +  (Xn+1) 

'  I  2a  2 


]IIP  „2  +  f2ifL  +  ^]lw  _wj! 

J  n  0  0  l  2Xna  2  y  n  0 


+  X0{8(llvn»  +  Iv  I)  +  2}IIPnv0  -  v  II  +  — «A(q  )v  -  A(q „)v  I*  . 

zee 

From  this,  the  uniform  bound  on  II v  L  lim  w  =  w.  in  H  and  assumptions 

n  n-.co  n  o  K 

(A)  and  (D)  allow  us  to  conclude  lim  v  =  v.  in  V  and  that  (3.7)  holds. 

n-*®  n  u 

An  analogous,  but  somewhat  simpler  argument  can  be  used  to  verify  (b). 
We  use  Theorem  2.2  to  show  that  for  n  e  Z  +  fixed,  lim  U  (t,s;q  )w 

—♦co  nv  ’  '^nv  m 

m 

Un(t,s;q0)wQ,  0  ^  s  U  S  T,  uniformly  in  t  for  t  €  [s,T]  whenever  wm,wQ  €  Hn 
with  lim  wm  =  wQ  in  H.  Clearly  lim  fn(-;qm)  =  f„(;q0)  in  L,(0,T;H  )  so  that 
we  need  only  to  show  that  for  some  XQ  >  0, 

Urn  J(X0;An(qm)  +  ul)w  =  J(\0;An(q0)  +  ul)w 

m 

in  H  whenever  lim  w  =  w  in  H.  Let  v_  =  J(X  *A  (q  )  +  ul)w  and  vn  = 

m_*co  ru  O  m  N  u*  n'^rrr  m  0 

J(X0;An(q0)  +  wl)w0-  Then  from  assumption  (B) 

X0al,Vm“V0H2  «  Vm~V0^  +  X0<A^n,)Vm  "  V0 - Vm- V0 > 

=  <(I  +  X0(An(qm)  +  wl))vm  -  (I  +  X0(An(q0)  +  •Jl))v0,v|B  -  vQ> 

"  lVm"Vol2  +  X0<A(%>V0  "  A^m)V0’Vm  "  V 
=  <Wm"W0’  Vm“V0>  ~  lVm_Vof2  +  “  A(<OV0’Vm  "  V 

S  Hwm-w0l  Jvm-v0ll  +  X0llA(q0)v0  -  A(qm)v0Uvm  -  v0l 


"vm  _voHr:  iwm  -  woi  +  “  ,A(%)vo  ~  A(qm)v0i.. 
xoa  a 


's, fOCrOZti.'Ov 


Assumption  (A)  and  iim  wm  =  wQ  in  H  yield  the  desired  result  and  the 
theorem  is  proved. 

Remark.  In  practice,  the  approximating  identification  problems  (IDn)  are  solved 
using  standard  iterative  search  techniques  (for  example,  steepest  descent,  Newton’s 
method,  etc.)  requiring  the  evaluation  of  $n(q)  for  q  €  Q  at  each  step.  This  in  turn 
requires  the  integration  of  the  finite  dimensional  initial  value  problem  (3.4),  (3.5). 
Once  a  basis  for  Hn  has  been  chosen,  the  solution  to  (3.4),  (3.5)  can  be  computed  using 
any  standard  numerical  integrator  for  ordinary  differential  systems.  Also,  the 
parameter  space  Q  and  the  admissible  parameter  set  Q  are  frequently  functional  in 
nature  and  infinite  dimensional.  When  this  is  the  case  the  set  Q  must  also  be 
discretized.  Suppose  that  for  each  m  =  1,2, ...,  I™:  Q  CQ  -*Q  is  a  continuous  map  with 
finite  dimensional  range  and  that  lim  I™(q)  =  q  with  the  convergence  uniform 

m-“ 

in  q  for  q  6Q.  Set  Qm  =  Im(Q)  (note  that  Q™  is  a  compact  subset  of  Q)  and 
consider  the  identification  problems  (ID™)  defined  to  be  the  problems  (IDn) 

with  Q  replaced  by  Qm.  It  is  clear  that  each  of  these  problems  admit  a 

solution  q^1  and  it  is  not  difficult  to  argue  that  there  exists  a  subsequence 

{qm* }  C  {q™}  with  lim  q™^  =  q,  q  a  solution  to  problem  (ID)  (see,  for  example, 
nk  n  j.k-*”  nk 

[4]).  Once  bases  for  Hn  and  the  range  of  I™  have  been  chosen,  problem  (ID™) 
involves  the  minimization  of  a  functional  over  a  compact  subset  of  Euclidean 
space  subject  to  finite  dimensional  constraints. 

Remark  (Nonautonomous  systems).  Theorems  2.1  and  2.2  remain  valid  for 
certain  classes  of  temporally  inhomogeneous  or  time  dependent  operators  A  = 
A(t).  To  be  more  precise,  the  family  of  operators  A(t):  X  -  2X  must  be 
m-accrctivc  on  X  for  almost  every  t  €  [0,T]  and  must  satisfy 
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(3.8)  |J(X;A(t))x  -  J(X;A(s))x|x  «  X|h(t)  -  h(s)|xL(|x|x) 

for  each  x  G  X,  every  X  satisfying  0  <  X  «  XQ  for  some  XQ  >  0,  some  h  G 
Lj(0,T;X),  some  continuous,  non  decreasing  function  L:  [0,°°)  -  (0,“)  and 
almost  every  t,s  G  [0,T]  (see  [8],  [9]).  (Note  that  for  simplicity  we  have  taken 
oj  =  0;  however,  the  discussion  to  follow  remains  valid  for  any  u  G  R.)  The 
primary  motivation  for  developing  the  framework  outlined  above  was  to 
define  readily  verifiable  conditions  on  the  operators  A(q):  V  -»  V*  that  if 
satisfied  would  (i)  also  automatically  be  satisfied  by  the  Galerkin 
approximation  An(q)  and  (ii)  lead  to  the  desired  convergence  of  solutions  to 
the  approximating  identification  problems  to  a  solution  to  problem  (ID).  The 
natural  assumption  to  add  to  (A)  -  (C)  that  certainly  satisfies  criterion  (i)  and 
that  could  conceivably  lead  to  an  estimate  of  the  form  (3.8)  in  H  is  that 

(3.9)  8A(t;q)v  -  A(s;q)vll.  «  |h(t)  -  h(s)|  L(|v|) 

for  each  v  6  V,  almost  every  s,t  G  [0,T]  and  some  h  G  1^(0, T;H)  and  some 
continuous  nondeercasing  L:  (0,®)  -  (0,“),  both  of  which  do  not  depend  upon 
q  G  Q.  Unfortunately,  however,  we  can  only  show  that  (3.9)  leads  to  an 
estimate  of  the  form 

(3.10)  |J(X;A0(t;q))u  -  J(  X;AQ(s;q))u|  S  /X  |h(t)  -  h(s)|  L(|u|) 

for  each  u  G  H.  Moreover,  it  is  not  clear  to  us  how,  or  if,  the  proof  of  the 
fundamental  Theorem  2.1  given  in  [9]  could  be  modified  so  that  (3.10)  would 
suffice.  We  have  explored  alternative  approaches  and  developed  other 
techniques  for  treating  the  nonautonomous  case  (for  example,  in  the  linear 
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4.  Applications  and  Examples 


We  briefly  describe  some  classes  of  systems  to  which  the  general 


framework  developed  in  the  previous  section  applies.  In  our  discussion  below 


we  consider  theoretical  aspects  only.  Implementation  questions  will  be  treated 


and  the  results  of  our  numerical  studies  will  be  reported  on  elsewhere. 


Example  4.1.  Linear  regularly  dissipative  operators.  The  approximation 


theory  for  inverse  problems  for  systems  involving  linear  regularly  dissipative 


operators  was  treated  in  detail  by  Banks  and  Ito  in,  and  is  the  central  focus 


of,  [2]  and  [3],  We  show  here  that  the  linear  theory  is  a  special  case  of  the 


nonlinear  theory  given  in  Section  3. 


Let  the  spaces  H,  V,  V*  and  Q  and  the  set  Q  be  as  they  have  been  defined 


above.  For  each  q  e  Q  let  a(q)(-,-)  be  a  sesquilinear  form  defined  on  V  x  V 


which  satisfies  the  conditions: 


(A1)  For  each  v  e  V  the  mapping  q  -  a(q)(-,v)  is  continuous  from  Q  C  Q 


into  V*.  That  is  given  €  >  0  there  exists  a  5  >  0  such  that 


sup  |a(qn)(u,v)  -  a(q)(u,v)|  <  e 

u€V  u 


whenever  d(qn,q)  <  6  where  d  denotes  the  metric  on  Q. 


(B1)  There  exist  an  w  €  R  and  an  a  >  0,  both  independent  of  q  €  Q,  for 


which  a(q)(v,v)  +  u|v|2  ?  oc II v II 2  for  every  v  €  V. 


(C1)  There  exists  a  constant  B  >  0,  independent  of  q  6  Q,  such  that 


|a(q)(u,v)|  S  BlluN  H v H  for  every  u,v  €  V. 


vvS> 


SSylSi 


When  conditions  (A1)  -  (C')  arc  satisfied  it  is  not  difficult  to  argue  that  for 
each  q  €  Q  an  operator  A(q)  6  £(V,V*)  can  be  defined  by 

[A(q)v](u)  =  <A(q)v,u  >  =  a(q)(u,v), 

u,v  €  V  and  that  A(q):  V  -  V*  satisfies  (A)  -  (C).  It  then  follows  from 
Theorems  2.1  and  3.1  that  there  exists  a  unique  nonlinear  evolution  system 
{UQ(t,s;q):  0  «  s  $  t  «  T}  on  H  corresponding  to  the  initial  value  problem 

u(t)  +  A0(q)u(t)  =  f(t;q),  0  <  t  S  T 
u(0)  =  u°(q) 

where  for  each  q  e  Q,  f(-;q)  €  Lj(0,T;H),  u°(q)  €  H  and  AQ(q):  Dom(A0(q))  C 
H  -•  H  is  the  restriction  of  A(q)  to  the  set  Dom(A0(q))  =  {v  e  V:  A(q)v  e  H}. 
The  operator  -A0(q)  is  the  infinitesimal  generator  of  an  analytic  semigroup 
(T0(t;q):  t  )  0)  on  H  (see  (18])  and  for  <£  e  H 

(4.1)  Uo(t,s;q)0  =  To(t~s;q)0  +  J  T0(t-T;q)f(T;q)dT. 

J  8 

It  can  be  shown  that  the  semigroup  {TQ(t;q):  t  ?  0)  admits  an  extension  (T(t;q): 
t  J  0}  which  is  an  analytic  semigroup  on  V*  with  generator  A(q):  V  C  V*  -*  V*. 
Also  the  restriction  of  (TQ(t;q):  t  2  0}  to  V,  call  it  (T(t;q):  t  £  0),  is  a  nalytic 
semigroup  on  V  with  generator  A(q):  Dom(A(q))  C  V  -»  V,  the  restriction  of 
A(q)  to  the  set  Dom(A(q))  =  (v  e  V:  A(q)v  e  V)  (see  [3],  [18]).  Consequently, 
with  appropriate  assumptions  on  f(  ;q),  the  evolution  system  (U0(t,s;q):  0  S  s 
S  t  «  T]  admits  an  extension  (U(t,s;q):  0  S  s  it  S  T)  which  is  an  evolution 
system  on  V*  and  a  restriction  (U(t,s;q):  0  ^  s  5t  S  T)  which  is  an  evolution 


system  on  V. 
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It  is  clear  from  (4.1)  that  when  A(q)  is  linear,  we  may  take  f(;q)  =  0  and 


consider  only  the  approximation  of  the  semigroup  {TQ(t;q):  t  £  0}.  For  each 


n  =  1,2,  ...  let  the  finite  dimensional  subspaccs  Hn  of  H  and  the  corresponding 


orthogonal  projections  P  be  as  they  were  defined  in  Section  3  and  assume 


that  condition  (D)  is  satisfied.  Denote  the  Galcrkin  approximations  to  A(q) 


(i.e.  the  restriction  of  A(q)  to  an  operator  from  Hn  into  H*  =  Hn)  by  An(q)  and 


set  T  (t;q)  =  exp(-tA  (q)),  t  ?  0.  Theorem  3.2  then  implies  that 


(4.2)  lim  |Tn(t;qn)Pnu°(qn)  -  T0(t;q0)u°(q0)|  =  0 


uniformly  in  t,  for  t  €  [0,T]  whenever  {qn}  C  Q  with  lim  qn  =  qQ  e  Q,  and 


the  mapping  q  -  uc(q)  is  continuous  from  Q  C  Q  into  H.  In  addition,  recalling 


that  we  required  that  Hn  C  V  for  all  n  =  1,2,  ....  an  inspection  of  the  proof  of 


Theorem  3.2  reveals  that  in  the  linear  case  with  the  existence  of  the  semigroup 


(T(f;q):  t  ?  0}  on  V,  we  may  apply  Theorem  2.2  with  X  =  V  and  conclude  that 


(4.3)  limlTn(t;qn)Pnu°(qn)  -  T(t;q0)u°(q0)l  =  0 


uniformly  in  t  for  t  €  [0,TJ  whenever  lim  qn  =  qQ,  u°(q)  €  V  and  the  map 

n~*°° 

q  -  u°(q)  is  continuous  from  Q  into  V  (see  also  [3]).  Then  for  (>6H,  setting 


Un(t,s;q)Pn<0  =  Tn(t-s;q)Pntf>  +  j  Tn(t-T;q)Pnf(T;q)dT 


under  appropriate  assumptions  on  f(;q),  (4.2)  and  (4.3)  continue  to  hold  with 


Tn(t;q),  T0(t;q),  and  T(t;q)  replaced  by  Un(t,s;q),  UQ(t,s;q),  and  U(t,s;q), 


respectively,  with  the  convergence  being  uniform  in  t,  for  t  €  [s,T],  Hence  the 


linear  theory  and  results  of  [3]  are  a  special  case  of  the  nonlinear  theory  of 


Section  3. 


if 
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Wc  note  that  in  the  context  of  the  identification  problem,  the  fact  that 

the  stronger  V-convcrgcnce  given  in  (4.3)  can  be  obtained  is  significant. 

Indeed,  (4.3)  permits  the  relaxation  of  the  continuity  assumption  on  the 

performance  index  4>  to  the  requirement  that  for  each  z  6  Z,  the  mapping  u  -* 

S’fu, z)  be  continuous  from  C([0,T];V)  into  R+.  This  can  have  the  effect  of 

significantly  enlarging  the  class  of  allowable  observations.  For  example,  in 

the  case  of  a  one  dimensional  parabolic  system  formulated  in  H  =  L2  with  V 

in  H1,  spatially  discrete  (i.e.  pointwise,  as  opposed  to  distributed  in  space) 

measurements  will  suffice  (see  [3]  and  [5]). 

Among  the  class  of  linear  regularly  dissipative  operators  which  arise  from 

a  form  satisfying  (A')  -  (C1)  arc  the  familiar  elliptic  partial  differential 

.  42+4+i 

operators  on  L_.  Briefly,  let  be  a  region  in  R  and  let  Q  =  x  L^ft).  Let 

2  n=l 

Q  be  a  compact  subset  of  Q  with  the  property  that  if  q=  {(a-,  L,  c):  i,j=l  €  Q, 
then  for  some  a  >  0  independent  of  q  6  Q. 
l 

I  a  (x)^  *  aK|2 
i ,  j  —  1  1J  1  J 


for  every  x  €  fi,  and  every  ^  €  R*.  For  q  e  Q  and  u,v  e  H1^)  set 


,  4 

(u,v)  =  {  I 

JnV)=i 


,  ,  3u(x)  9v(x) 

Vx>  ax.  — 

l  J 


+  I  b(x)  v(x)  +  c(x)u(x)v(x)]dx 

i=t  1  9x.  J 


with  H  =  L2(fl)  and  V  any  closed  subspace  of  H*(n)  containing  H*(fi),  it  can  be 
shown  (see  [18])  that  a(q)(-,-)  satisfies  (A’)  -  (C’).  The  operator  A(q)  is  given 


When  30  is  sufficiently  smooth,  A(q)  is  the  elliptic  operator  given  by  (4.4),  and 
V  is  chosen  to  be  either  Hq(0)  or  H1(0),  the  equation  (3.1)  becomes  a  parabolic 
partial  differential  equation  with  either  Dirichlet  or  Neumann  boundary 
conditions. 

For  H  =  L2(0)  and  V  a  subspace  of  H1(0),  choosing  the  approximating 
subspaces  to  be  the  span  of  an  appropriate  collection  of  first  order  spline 
functions  will  typically  satisfy  assumption  (D)  (see  [15]  and  Example  4.2 
below). 

Example  4.2.  Nonlinear  Elliptic  Operators.  Let  0  be  a  bounded  region  in  R® 
with  smooth  boundary  r  =  30.  For  a  =  (o1,  ...,  oc^)  a  multi-index,  let  |a|  =  ot  + 
a2  +  •  -  +  and  denote  the  orth  order  generalized,  or  distributional 

derivative  of  a  function  u  by  D®u;  that  is, 

a,  a. 

„  3  1  3  ® 

Dau(x)  =  — —  •  •  -  — —  u(x),  x  €  0  . 

3x“®  3x“® 

Let  m  be  a  nonnegative  integer  and  let  8u  denote  the  vector  valued  function 
of  length  N  =  (®Jm)  whose  components  are  all  of  the  partial  derivatives  of  u 
of  order  greater  than  or  equal  to  zero  and  less  than  or  equal  to  m. 

For  each  multi-index  a  with  |aj  $  m,  let  (x,£)  -*  aa(x,0  be  a  real  valued 
function  defined  on  0  x  RN  which  is  measurable  in  x  and  continuous  in  We 
assume  that 

(I)  there  exist  a  g  €  L2(fi)  and  a  positive  constant  7  such  that 


(4.5) 


|aa(x,;)|  «  7(KI  +  g(x)) 

for  almost  every  x  €  0,  each  t;  G  RN  and  all  ot  with  jaj  ?  m,  and 
(2)  there  exists  a  positive  constant  X  such  that 


(4.6) 


E 

Jccj 


(aa(x,0  -  aa(x,n))(;a  -  na)  $ 


x  E 

|a|Sm 


-  V2 


for  almost  every  x  G  Q  and  all  t;,n  G  RN. 

Let  H  =  L2(n)  and  let  V  be  any  closed  subspace  of  Hm(fi)  which 
contains  H^(fl).  Define  the  operator  A:  V  -  V*  by 

(4.7)  (Au)(v)  =  E  [  a  (x,6u(x))Dav(x)dx, 


for  u,v  G  V.  The  operator  A  given  by  (4.7)  is  the  distributional  form  of  the 
formal  differential  operator 


(4.8)  (Au)(x)  =  E  (-l^D0^  (x,6u(x)). 

|oc]Sm 

A  differential  operator  of  the  form  (4.8)  is  referred  to  as  a  nonlinear  elliptic 
operator  and  the  partial  differential  equation 

(4.9)  ~  (t,x)  +  E  (-l)aDaa  (x,5u(t,x))  =  f(t,x) 

dt  |a|Sm 

is  said  to  be  of  nonlinear  parabolic  type.  When  V  =  H™(ft),  a  solution  in  V* 
to  the  abstract  equation 


u(t)  +  Au(t)  =  f(t) 

with  A  given  by  (4.7)  corresponds  to  a  variational  solution  to  (4.9)  which 
satisfies  Dirichlct  boundary  conditions.  When  V  =  Hm(fl),  a  variational 


solution  to  the  Neumann  problem  is  obtained.  Note  that  in  the  linear  case  we 
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aa(x,6u(x))  =  E  aa'6(x)D6u(x). 

Ini  $ m 


Under  the  assumptions  above,  it  is  not  difficult  to  show  that  A  given  by 


(4.7)  is  hemicontinuous  and  satisfies  conditions  (B)  and  (C)  given  in  Section  3. 


With  an  appropriate  choice  of  the  space  Q  and  the  set  Q,  condition  (A)  can  be 


satisfied  as  well. 


A  quasilinear  model  for  heat  conduction  or  mass  transfer  in  which  the 


heat  or  mass  flux  is  a  function  of  the  temperature  or  mass  fraction  gradient 


discussed  in  [16]  and  [17]  leads  to  a  nonlinear  elliptic  operator  and  a  nonlinear 


parabolic  partial  differential  equation  of  the  forms  (4.8)  and  (4.9),  respectively 


with  m  =  1.  Let  ft  be  a  bounded  region  in  R*  with  smooth  boundary  and  let 


Q  =  LJO  x  R*).  Let  Q  be  a  compact  subset  of  Q  with  the  property  that 


q  <=  Q  if  and  only  if  the  mapping  Q  -  q(x,0  is  C1  on  R*  for  almost  every 


x  €  O  and  there  exists  a  X  >  0  (which  does  not  depend  on  q)  such  that 


(4.10)  e^qOc.Ol^g-U  ~  »»)  +  q(x,9)U.  -  n.)  *  X({.  -  »»,), 


for  i  =  1,2,  ...,  *,  almost  every  x  €  O  and  all  6,5, n  €  R*.  (When  1  =  1,  the 


function  q(x,0  =  q(()  =  (1  -  ,5e'*  )  satisfies  (4.10).) 


Let  H  =  L,(0)  and  let  V  be  any  closed  subspace  of  H^fl)  which  contains 


Hg(fl).  Then  V  c  H  C  V*  and  for  each  q  e  Q  define  A(q):  V  -  V*  by 


(4.11)  ( A(q) 


u)(v)  =  q(x. 


Vu(x))Vu(x)-  Vv(x)dx 


for  u,v  €  V.  Note  that  for  each  q  €  Q  the  operator  given  by  (4.11)  is  of  the 


form  (4.7)  with 
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(4-12)  aa(x,6u(x))  =  q(x,Vu(x))Dau(x) 


for  x  €  Q  and  all  a  with  |aj  =  1  and  aa  =  0  for  |aj  =  0.  The  nonlinear 


parabolic  partial  differential  equation  (4.9)  takes  the  form 


—  (t,x)  -  V-q(x,vu(t,x))Vu(t,x)  =  f(t,x),  t  >  0,  xen. 

ot 


Taking  II  •  II  to  be  the  usual  norm  on  H^ft),  it  follows  that 


II A(q0)u  -  A(qx)u II*  S  |qQ  -  q^  Hull 


for  each  u  €  V  and  q^qj  €  Q.  Since  Q  is  a  compact  subset  of  x  R*),  it  is 


easily  verified  that  aa  given  by  (4.12)  satisfies  a  growth  condition  of  the  form 


(4.5)  with  7  and  g  independent  of  q  €  Q.  An  application  of  the  mean  value 


theorem  together  with  assumption  (4.10)  imply  the  existence  of  a  X  >  0, 


independent  of  q  €Q,  for  which  (4.6)  holds.  Consequently  the  conditions  (A),  (B), 


and  (C)  given  in  Section  3  are  satisfied,  and  our  general  theory  can  be  applied. 


With  regard  to  approximation,  polynomial  spline  function  based  Galerkin 


subspaces  can  often  be  shown  to  satisfy  condition  (D).  For  example,  when 


S  =  1  and  O  =  (0,1)  in  the  nonlinear  heat  conduction/mass  transfer  example 


discussed  above,  the  subspaces  FIn  can  be  chosen  as  the  span  of  the  linear 


B-spline  ("hat")  functions  with  respect  to  the  uniform  mesh  {0,  1/n,  2/n,  ...,  1} 


appropriately  modified  to  satisfy  stable,  or  geometric,  boundary  conditions. 


Familiar  error  estimates  for  interpolation  and  the  Schmidt  inequality  can  then 


be  used  to  verify  that  condition  (D)  is  satisfied  (see  [5]).  Generalization  to 


higher  dimensions  is  possible,  and  can  often  be  achieved  via  tensor 


products  of  one  dimensional  elements  (see  [15]). 
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5.  Concluding  Remarks 


Wc  have  developed  a  general  abstract  approximation  framework  for  the 
identification  of  nonlinear  distributed  parameter  evolution  systems.  The  class 
of  systems  to  which  our  theory  applies  are  those  whose  dynamics  can  be 
described  by  a  nonlinear  operator  which  satisfies  conditions  that  arc  the 
natural  nonlinear  extensions,  or  analogs,  of  the  properties  of  regularly 
dissipative,  or  abstract  parabolic,  linear  operators.  The  approach  we  have 
taken  is  based  upon  the  defining  of  a  sequence  of  approximating  finite 
dimensional  identification  problems  in  which  the  systems  to  be  identified  are 
Galerkin  approximations  to  the  original,  underlying,  infinite  dimensional 
nonlinear  dynamics.  Under  a  weak  continuity  assumption  with  respect  to  the 
unknown  parameters  to  be  identified,  equi-boundedness  and  equi- 
monotonicity  conditions,  and  an  approximation  assumption  on  the  Galerkin 
subspaces  (all  of  which  are  readily  verified  for  wide  classes  of  nonlinear 
distributed  systems  and  finite  element  subspaces),  we  are  able  to  demonstrate 
that  solutions  to  the  approximating  problems  exist,  and,  in  some  sense, 
approximate  (i.e.  subsequential  convergence)  solutions  to  the  original  infinite 
dimensional  identification  problem.  We  have  shown  that  the  linear  theory 
presented  in  [2]  and  [3]  is  a  special  case  of  our  nonlinear  framework  and  that 
our  results  are  applicable  to  a  reasonably  wide  class  of  nonlinear  elliptic 
operators  and  corresponding  nonlinear  parabolic  partial  differential  equations. 
In  particular,  we  have  considered  application  of  our  theoretical  framework  to 
a  quasi-linear  model  for  heat  conduction  or  mass  transport. 
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Thc  general  approximation  result  for  nonlinear  evolution  systems 
discussed  in  Section  2  is  applicable  to  a  much  broader  class  of  nonlinear 
dynamical  systems  than  we  subsequently  treated  in  Section  3.  For  example, 
this  class  of  systems  would  include  those  with  dynamics  described  by  set 
valued  maps  or  multifunctions,  and  (after  minor  modification  to  the 
general  theory)  time  dependent  or  nonautonomous  operators.  We  are 
currently  investigating  these  features  of  the  general  approximation  theory 
in  the  context  of  parameter  estimation  problems.  Also,  we  would  like  to 
be  able  to  weaken  the  somewhat  restrictive  strong  monotonicity  condition. 
Any  progress  that  we  might  make  in  these  efforts  would  have  the  potential 
to  significantly  enlarge  the  class  of  nonlinear  systems  to  which  our  theory 
and  framework  would  apply.  Finally,  extensive  numerical  or 
computational  studies  designed  to  demonstrate  the  feasibility  and  point  out 
the  limitations  of  our  schemes  and  general  approach  are  currently 
underway  and  will  be  reported  on  in  a  forthcoming  paper. 
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